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We investigate how thin sheets of arbitrary shapes morph under the isotropic in-plane expansion 
of their top surface, which may represent several stimuli such as nonuniform heating, local swelling 
and differential growth. Inspired by geometry, an analytical model is presented that rationalizes how 
the shape of the disk influences morphing, from the initial spherical bending to the final isometric 
limit. We introduce a new measure of slenderness 7 that describes a sheet in terms of both thickness 
and plate shape. We find that the mean curvature of the isometric state is three fourth’s the natural 
curvature, which we verify by numerics and experiments. We finally investigate the emergence of a 
preferred direction of bending in the isometric state, guided by numerical analyses. The scalability 
of our model suggests that it is suitable to describe the morphing of sheets spanning several orders 
of magnitude. 
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Cylindrically curved thiu structures result from the 
uauo-scale fabricatiou of semicouductor uauotubes HHa, 
aud the uouuuiform heatiug [5], local swelliug [6], aud 
differeutial growth [ 7 ] of thiu sheets. For lamiuated com¬ 
posites [8], electrolytic thiu film depositiou aud 

coucrete slabs m , this cyliudrical curliug preseuts au eu- 
giueeriug challeuge, while receut work has utilized it as a 
mechauism for stimuli respousive self-assembly nans]. 
The leugth scales of these examples rauge from the 
uauometer to the meter, suggestiug that geometry dom- 
iuates the deformatiou processes. Mechauically, these 
structures are bilayer disks iu which oue layer isotrop¬ 
ically expauds relative to the other. Iu this Letter, we 
show that iu the asymptotic limit of large bilayer growth, 
auy arbitrarily shaped disk will adopt a cyliudrical shape 
whose meau curvature is three-fourth’s the uatural cur¬ 
vature. We preseut au aualytical model that captures 
both the bifurcatiou from spherical to cyliudrical curva¬ 
ture [14] aud the isometric limit, verified by uumerics aud 
experimeuts. 

Let us cousider a thiu body subjected to a dome-like 
(elliptic) uatural curvature - this may be thought of as 
the result of swelliug, growth, or heatiug of the top sur¬ 
face. Withiu the coutext of uou-Euclideau plates HSHH], 
the body may be modeled as a shell haviug the followiug 
first aud secoud uatural fuudameutal forms 


a = A 


2 

o 




( 1 ) 


The uatural stretch Aq aud uatural curvature Kq repre- 
seut the lateral distauces aud curvatures that would make 
the sheet locally stress-free. However, as every surface 
iu space must fulfill the Gauss-Codazzi-Maiuardi equa- 
tious, it is uot usually possible for the sheet to realize 
the uatural forms. Gouversely, if a beam is cut from the 
sheet, its iuhereut ID geometry will allow it to adopt 
a three-dimeusioual shape with lougitudiual axis stretch 
aud curvature equal to Aq aud Kq, respectively, without 



FIG. 1. (a) A bilayer disk and the beam cut from it. (b) 
Isometric limit, (c) Measuring the natural curvature of the 
disk as the realized curvature of the beam, (d) Isometric 
states obtained by residual swelling of disks having different 
shapes. 


any need of satisfying additional constraints. This strong 
connection between the sheet and the beam provides a 
simple means to evaluate the natural forms of the for¬ 
mer by measuring the deformed shape of the latter [see 
figure [^(c)]. 

To fully exploit this idea, we fabricated bilayer disks 
by spin-coating the two layers made of polyvinylsilox- 
ane (PVS) Zhermack Elite Double 32 and 8, respectively, 
as shown in figure (a-b). Eree chains flow from the 
softer (pink) to the stiffer (green) polymer by residual 
swelling m, corresponding to a dome-like natural cur¬ 
vature, and a conformal natural stretch. Mechanically, 
this corresponds to a bilayer disk where the top layer 
has been stretched homothetically by a factor and 
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then bound to the bottom layer. This results in a dis¬ 
continuous three-dimensional natural metric g that may 
be approximated as g = a — 2zh, where 2 ; is the coordi¬ 
nate along the thickness. Approximate expressions for a 
and b as functions of pre-stretches were derived in [ 18 ], 
albeit the natural curvature was then measured experi¬ 
mentally as the curvature of very small strips cut along 
principal directions. We overcome the lack of a predictive 
expression for the natural curvature by imaginarily cut¬ 
ting a beam from the sheet, and determining its deformed 
shape analytically by following m as Ao = Ao(A,m,n) 
and Ko = hZo{X,m,n)/h^ where m and n are top-to- 
bottom thicknesses and Young moduli ratios, respec¬ 
tively, and h is the thickness of the sheet. Therefore, 
the natural forms are determined analytically as func¬ 
tions of A and the material and geometrical ratios. 

A comparison of the energies for stretching and bend¬ 
ing suggests that a sheet exposed to such stimuli should 
adopt an isometric deformation in the limit of large 
stretch. We analytically describe the realization of this 
asymptotic isometry^ corresponding to a zero Gaussian 
curvature K and a nearly stretch-free sheet. The dimen¬ 
sionless energy of the shell may be written as m 

U= y*[(l — z/)|a — ap + z/tr ^(a — a)]y^(iA+ 

+ y J[{l-I^)\h-h\‘^ + ptr‘^{h-h)]^/\i\dA, 

( 2 ) 

where a and b are the first and second fundamental forms 
of the mid-surface of the deformed sheet, respectively. In 
the isometric limits we have a = a and the stretching 
energy is zero. To find the deformation of the disk we 
then minimize the bending energy under the constraint 
that the mid-surface be fiat, since a = a = A^I. We 
write b in cartesian coordinates (iz, v) so that the second 
fundamental form is Ldv? ^2Mdudv^Ndv‘^ ^ and impose 
the constraint of a flat mid-surface by enforcing that its 
Gaussian curvature A~^{LN — M^) be zero through a 
Lagrange multiplier. 

For the characterization of the isometric limits we as¬ 
sume that the energy density is homogenous. Gonse- 
quently, for the total energy to be minimized, it suffices 
to minimize its density augmented by the constraint on 
the Gaussian curvature imposed by the Lagrange mul¬ 
tiplier. The minimization gives L N = Ko{ 1 -\- u) 
and LC = 0, a solution solely in terms of the principal 
invariants of the curvature tensor, and independent of 
the shape of the sheet: if a homogeneous energy density 
is assumed, no information on the bending direction can 
be obtained [20| . The result of the minimization may be 
rewritten in terms of the mean curvature H of the sheet 
as iL = = |A-2(L + AT). If the 


material is incompressible (z^ = 1 / 2 ), it becomes 


that is, in the limit of large stretch, the mean curva¬ 
ture of the sheet is three-fourth of the rescaled natural 
curvature z^J = z^o/Aq, regardless of the shape of the 
sheet. Equivalently, the sheet morphs into a cylinder of 
radius R = 1/{2H). Since Kq ^ 1/h whereas Aq does not 
depend on the thickness, we derive H ^ 1/h. 

To verify our analytical prediction, we carried out 
experiments and numerical simulations. In the exper¬ 
iments, we measured the curvatures of disks of sev¬ 
eral shapes and the corresponding beams during resid¬ 
ual swelling [figure (a-b)] by image analysis performed 
using Matlab. Numerical simulations were performed to 
solve the geometrical problem within the context of fi¬ 
nite (incompatible) tridimensional elasticity with large 
distortions using a Neo-Hookean incompressible mate¬ 
rial model [21] implemented in the commercial software 
GOMSOL Multiphysics. The sheets were made of two 
layers: as the top one was pre-stretched by A“^, it was 
subjected to a distortion field = A(ei G ei -h 62 ( 8 ) 62 ) + 




FIG. 2. (a) Dimensionless mean curvature versus dimension¬ 
less rescaled natural curvature: the three-fourth law excel¬ 
lently captures numerical and experimental data regardless of 
initial shape. Numerical data (blue) refers to circle, square, 
leaf-like, triangle, cross and ellipse, each for different thick¬ 
nesses, and experimental data (red) refers to a square, cross, 
triangle, and circle, (b) Dimensionless mean curvature versus 
A: the isometric limit is expressed in terms of the stretch of 
the top layer. Three different master curves are shown, each 
corresponding to different thicknesses ratios m. (c) Isometric 
shapes from COMSOL with the color code representing the 
mean curvature. 
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es ( 8 ) es, whereas the bottom one was not pre-stretched 
and therefore subjected to = I. We simulated several 
shapes and thicknesses, and verified that the mean cur¬ 
vature H was homogeneous apart from boundary layers 
near the edges. Then we plotted the product as a 
function of the dimensionless natural curvature or 
the stretch A via the ID analytical model, in the limit 
of large stretch. The numerical data referring to differ¬ 
ent thicknesses and the experimental data collapse to a 
single master curve thus verifying the analytical predic¬ 
tion that H 1/h and the three-fourth law excellently 
predicts the mean curvature of the sheet regardless of its 
shape as shown in figure(a). Figure (b) shows how 
the dimensionless mean curvature varies with the stretch 
of the top layer and how it is affected by the thicknesses 
ratio m, providing a design rule for bilayer sheets; all the 
data corresponding to different thicknesses and shapes 
collapses to master curves. 

To predict and explain the transition that turns a 
spherical growth into a flat state [ 20 ], we start from the 
energy in equation (§ and compare the energies of the 
spherical and isometric states, without any a priori as¬ 
sumptions of a stress state or displacement. These as¬ 
sumptions usually rely on the consideration of an Airy 
stress function, which is known only for a very small num¬ 
ber of sheet’s shapes, and on the account of a linearized 
version of the Gauss’s theorem nni EH ng. Away from 
large stretches, the two principal curvatures are equal to 
each other and homogeneous throughout the sheet. How¬ 
ever, the morphing into a spherical cap {L = N) becomes 
too costly for the sheet when the stretch, or the natural 
curvature, reaches a critical value. While for small val¬ 
ues of /^o, it is convenient for the sheet to stretch and 
increase its Gaussian curvature, above a critical thresh¬ 
old of /^o, it is cheaper to morph isometrically. Assuming 
a metric with constant Gaussian curvature in Gaussian 
normal coordinates as in HZ], the dimensionless stretch¬ 
ing energy may be written as Us = j r^dA, 

where r is the undeformed radial coordinate, and we 
define S = (2/9 f r^dA/A)^!^ as the shape factor [24] . 
Also, as L = A", the dimensionless bending energy 
is U\y = h‘^A~^A(L — where A is the undeformed 

area of the incompressible sheet. As A~^L^ before 

bifurcation, the total dimensionless energy for the spheri¬ 
cal cap is = (1/9)T^A“^ J r^dA + h?AA~‘^(L — 
while the total energy of the isometric shape is equal 
to = (l/4)h^AA“^b:Q, which becomes smaller than 
the energy of the spherical cap above a critical value of 
the natural curvature m- This analytical reasoning sug¬ 
gests that a bifurcation between the two principal curva¬ 
tures must occur at a critical z^ob to drive the sheet into 
an isometric shape. 

As our model does not rely on any a priori assumption 
on the stress state or displacements, it is suitable to study 
the pre-bifurcation morphing of the sheets, regardless of 
their shapes. Minimization of the total energy U^b with 



FIG. 3. (a) Dimensionless natural curvature at bifurcation 
versus 7 : the analytical prediction given by (solid line), 
in agreement with the numerical solution (dashed) of 
perfectly predicts the finite element simulation results (sym¬ 
bols) obtained for different shapes (identified by the shape of 
the symbol) and thicknesses, (b) Pre-bifurcation behavior in 
terms of L('Ko) for different shapes (circle, square, rectangle) 
obtained by finite element simulations (symbols) and analyt¬ 
ics (solid lines). The two straight lines identify the upper and 
lower bounds corresponding to a bending-dominated solution 
and the bifurcation point, respectively, (c) Deformed shapes 
of a thin circular disk from COMSOL. 

respect to L yields 

P + 74(I-Ko) = 0 , (4) 

where curvatures have been nondimensionalized by 1 /h 
and denoted with bars. The geometrical dimensionless 
small parameter ^ = h/S encloses the shape of the sheet 
in the equilibrium equation and may be computed easily 
by its definition. The smaller 7 is, the more slender the 
structure is in terms of both thickness and plate shape. 
Equation 0 can be solved analytically and it admits 
just one real solution L(ko) or, equivalently via the ID 
model, L{X). 

When the natural curvature reaches a value such that 
the energy of the spherical cap coincides with the energy 
of the isometric state, bifurcation occurs; the equality 
between the two energies may be written as 

2 “ 2Lb7^ob + 4 ^ob) = 0 5 (5) 

where the subscript b denotes quantities at bifurcation. 
This is a perturbation problem in the smallness pa¬ 
rameter 7 that multiplies the bending factor of equa¬ 
tion Seeking a perturbation expansion solution in the 
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form I /5 = 1 / 50 + 7^^61 + 0 ( 7 ^), we obtain Lb — /^ob/ 2 , at 
the leading order. Moreover, by Taylor expanding equa¬ 
tion in 7 up to 0 ( 7 ^^/^), we derive a shape-dependent 
formula for the bifurcation natural curvature 

J20 + UV2 2 

Differently from other works in the literature, this for¬ 
mula remarkably predicts the natural curvature at bi¬ 
furcation for thin sheets having different shapes. While 
curvatures scale as 1 /h as we showed, this formula pre¬ 
dicts that the bifurcation natural curvature scales as h. 
The analytical solution evaluated at bifurcation 

through equation §, yields Li^ = 0.63Kobi analyti¬ 
cal relation that does not depend on 7 , and refines our 
previous asymptotic estimate. Figure (a) shows the 
dimensionless natural curvature at bifurcation for dif¬ 
ferent geometries (different 7 ): the analytical model ex¬ 
cellently predicts the threshold for different shapes and 
thicknesses that have been simulated in COMSOL pro¬ 
vided that the shape is convex. The dashed blue curve 
corresponds to the numerical solution of equation (§, 
which justifies the asymptotic analysis that led to equa¬ 
tion (§. The model is also capable of distinguishing the 
pre-bifurcation behavior that sheets of different convex 
shapes exhibit as figure (b) shows. For clarity, we 
reported the pre-bifurcation solution for three shapes. 
As 7 ^ 00 , bending dominates stretching and equa¬ 
tion Q yields L = Kq. This is also an upper-bound 
for all the shapes as L = 7 ^/ ( 7 ^ + 3L ) that is maxi¬ 
mum at Lo = 0 where it is equal to 1 regardless of the 
shape: the realized curvature never exceeds the natural 
curvature. Also, we plotted the relation between real¬ 
ized and natural curvatures at bifurcation, which sets a 
lower-bound for L. Figure [^(c) shows the numerical de¬ 
formed shape of a circular disk as the stretch increases: 
the disk morphs from a flat state to a spherical cap be¬ 
fore bifurcating into a cylinder. The analytical model is 
not capable of describing the spherical bending and the 
bifurcation of concave shapes, whose topological defects 
presumably affect their behavior [26]. Notably, the three- 
fourth law holds regardless of the shape, whether convex 
or concave. 

Finally, we investigate the emergence of a preferred di¬ 
rection of bending following bifurcation by running sev¬ 
eral simulations with different meshes. While it was 
shown that the principal curvatures of an ellipse have no 
preferred orientation [ 20 |, here we show that nearly all 
shapes select a particular direction for their generators 
immediately after symmetry is broken. Cylindrical gen¬ 
erators orient themselves along an axis which preserves 
reflection symmetry. If multiple planes of reflection sym¬ 
metry exist for a particular shape, bending will occur 
along the longest length of the disk (Fig. [^). For in¬ 
stance, a square of length a and width b has two unique 



FIG. 4. (a) Numerical statistics of preferred generators (blue) 
for different shapes, (b) Logistic analysis of the orientation 
data of the transition from a square to a rectangle. 


planes of reflection symmetry, and has been observed to 
bend along its diagonal m- We show with a series of 
simulations for a = 0.02 m and h = 0.4 mm that as a /6 is 
increased, and the disk transitions from square to rectan¬ 
gular, the orientation of the generators transitions from 
the diagonal to the short length 6 . Curiously, this transi¬ 
tion does not occur immediately once a /6 > 1 , as a logis¬ 
tic analysis of the orientation data identifies a transition 
at a /6 = 1.028 ± 0.0025. A preferred bending direction 
in rectangular sheets has been attributed to boundary 
layers effects [28| that lower the energy in the disk m- 
Our results help generalize these findings for any arbi¬ 
trary shape. The notion of a preferred bending direction 
upon bifurcation does not contradict Mansfield’s classi¬ 
cal results, as the cylindrical generators can be rotated to 
any arbitrary angle while remaining in equilibrium when 
perturbed via mechanical or chemical stimuli [ 6 ]. 

In summary, we have derived a theoretical analytical 
model that describes the morphing of disks subjected to 
dome-like natural curvatures, regardless of their shape. 
We also showed the emergence of an isometric limit with 
analytics, numerics and experiments, underlying the im¬ 
portance of isometries in mechanics m- The capability 
to morph 2D shapes into 3D shells by in-plane and 
transverse residual swelling opens intriguing avenues to¬ 
wards the precise design of soft structures. 
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